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Fermat's little theorem states that for p prime and a G Z, p divides — a. 
This result is of huge importance in elementary and algebraic number theory. 
For instance, with its help we obtain the so-called Frobenius automorphism of 
a finite field Fpn over ¥p. 

This theorem has many interesting and sometimes unexpected proofs. One 
classical proof is based upon properties of binomial coefficients. In fact, {d + 



l)P~dP~l^ Er=i {'^d\ Since (^) = is divisible by p for 1 < i < p - 1, 



then {d + 1)p — — 1 is divisible by p. Summing this over d = 1,2, a — 1, 
we obtain the desired result. Another classical proof is based upon Lagrange's 
theorem, which states that the order of an element of a finite group divides the 
group order. Applying this theorem to the multiplicative group of a finite field 
Fp we obtain the result immediately. Several other proofs can be found at p]. 
Nevertheless, in all of these proofs one or another analogue of the Euclidean 
algorithm (hence arithmetic) is being used. 

In this short note we present a curious proof which was found as a side result 
of another, unrelated problem (which is the case, maybe, with many such "cu- 
rious" proofs). Surprisingly, arithmetic, algebra, and the properties of binomial 
coefficients do not manifest at all. 

Let f{x) — 1 ^ X ~ dx^ + X]fc>3 (^kx'^ be any formal power series in Q, with 
coefficients in Z. It is well known that this series can be represented in a unique 
way as a formal product of the following form: 



where the coefficients nik are integers. This result can be found in 1 , but the 
proof is simple and straightforward. In fact, for k = 1 and k = 2 we have a 
unique choice mi = 1 and m2 = d. Suppose N > 3 and we have already chosen 



ruk for k < N ~ 1. Then n^i^(l ~ rnkx'^) = 1 — x — dx^ + J2k=3 '^^x^ + 



Cx + "higher terms" , where C is a certain integer which depends only on nik 
for 1 < A: < — 1. Therefore, the unique choice for tojv is rriN = C — qn- In 
a similar fashion, since = 1 + x + {d + l)x'^ + J2k>3 ^kX^ is also a formal 
integer power series, it can be represented in a unique way as a product 



f{x)^\[{l~mkx''), 



1 



(l + a;)(l + (d+l)a;2)[|( 1 — UkX ) 



fe>3 



where nk are integers as well, ni = —1, and n2 = —{d + 1). 

Recall that the logarithmic derivative of a power scries g{x), denoted by 
{\ng{x)y, is defined to be the power series g'{x)/g{x). It is not hard to prove 
that for any two formal power series g{x) and h{x), (hig{x)-h{x)y = (In g{x))' + 
i\nh(x))' . Indeed, this property reduces to the Leibniz rule 

{g{x)-h{x))' =g'{x)h{x)+g{x)h'(x). 

This is verified simply by comparing the corresponding coefficients. Note also 
that the binomial theorem is not used in the proof. 

Now take the formal logarithmic derivative of f{x). We obtain: 

kmkx^ AfV^ , N 



In a similar fashion, 



■' ^ ' N>1 s\N 

Therefore, we have interesting identities among the terms of two infinite se- 
quences: 



N ^ . N 



^^N/s—^-2^T^N/s—, N en. (1) 

|Af s\N 



We can easily prove by induction that this implies mfc = — n.^. for odd k, but 
not for the terms with even indices! Thus, a consequence of this reasoning is 
the fact that any infinite sequence of integers {mk,k € N} with mi = ±1 has 
an "inverse" sequence of integers {nk,k G N} with m = ^1. Consequently, all 
such sequences split into mutually inverse pairs. It is rather tempting to try to 
express an inverse of a certain sequence for which the infinite product has a rich 
mathematical content. For example, let us take nik = 1 for fc e N. Hence, we 
have a product 



{x,x)oc = 



fc=l 



It is well known that (x, x)^ = ^ p{n)x"' , where is Ramanujan's partition 

n=0 

function. Using the recurrence ([T]) we can compute the sequence rik — —rik- As 
mentioned, Uk — 1 for k odd, and terms of this sequence with even indices begin 
with 

2, 4, 0, 14, -4, -8, -16, 196, -54, -92, -184, 144, -628, -1040, -2160, 41102... 
Therefore, 



Ep(")a;" = Y[{l + nkx''' 



n=0 k=l 
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Let us return to our case. Recall that m2 = d and 712 = -~{d+ 1). Hence, 
when N — 2p, where p > 2 is a prime, ([T]) reads as: 

2p ■ map +p-ml + 2dP + 1 = -2p ■ n2p - p ■ nl + 2{d +1)^-1. 

Thus, p divides {d+l)P — dP — 1. Summing this over d = 1,2, ...,a — l, we finally 
obtain p\aP — a. Quite unexpected! 

Likewise, expand the following function into a formal infinite product: 

00 00 

fix) ^i-x-j2 ^^"2:"+! = n (1 - «"^")- 

n— 1 n— 1 

Since f{x) — — j^z^^^j after taking the logarithmic derivative, we obtain: 

-x(in/(x)) ^ y: e-'^e^at- 

N=l N>1 s\N 

As a direct consequence, Op — '^'''^^^ ^'^ ~^ , which implies that '^'''^^^ ~^ is 
an integer. Possible variations on this theme unexpectedly produce other con- 
gruences and identities. Recall that a prime number p is said to be a Wieferich 
prime if and only if 2^^^ = 1 (mod p^). Examples are p = 1093 and p = 3511, 
with no others in the range p < A ■ 10^^. In the last example with c? = 1, all 
the numbers Op = ^ ~^ appear simultaneously in the infinite product defining 
i^^, and as the proof of the algorithm used to expand a formal power se- 
ries into an infinite product suggests, strangely enough, the coefficients ajv are 
defined inductively on N without a distinction between prime and composite 
values of N . Possibly, more profound research of this product could clarify our 
understanding of these exceptional Wieferich primes. 
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